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We investigate the correlation between the neutrino mass limit and dark energy with time evolving 
equation of state. Parameterizing dark energy as w = wo + wi * z/(l + z), we make a global fit 
using Markov Chain Monte Carlo technique to determine wq,Wi, neutrino mass as well as other 
cosmological parameters simultaneously. We pay particular attention to the correlation between 
neutrino mass Em,, and Wi using current cosmological observations as well as the future simulated 
datasets such as PLANCK, SNAP and LAMOST. 
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I. INTRODUCTION 

o 

With the accumulation of high quality observational data of cosmic microwave background (CMB) [HQ, large scale 
structure (LSS) of galaxies [1,13, [B[ and Supernovae Type la (SN la) [1,0,11], we have been pacing into a new era of 
£>Y precision cosmolo gy. In 1998, the analysis of the luminosity-redshift relation of SN la revealed the current acceleration 
of our universe [9l llOfl. This acceleration attributes to a new form of energy, dubbed dark energy (DE), with negative 
pressure and almost not clustering whose nature remains unveiled. The simplest candidate of dark energy is the 
cosmological constant (CC) however it suffers from the fine tunning and coincidence problem s [T lL[T3|. To ameliorate 
these dilemmas some dynamical dark energy models such as quintessence jl3j |. phantom [l4|. k-essence fl5| and 
I quintom whose equation of state (EOS) can cross -1 during evolutio n II 611 h ave been proposed and studied thoroughly 
both theoretically and phenomenologically in the literature [H, G3, EmI, H HI SI Si H SI, SI, S3, SI, SH . 
Neutrino physics, especially the mass limit, is another challenge of modern science. The combined analysis of all 
. currently available data of neutrino oscillation experiments implies the crucial mass differences in the neutrino mass 
hierarchy, say, Am 2 2 ~ 7 x 10 _5 eV 2 , Am| 3 ~ 2.6 x 10~ 3 eV 2 for solar neutrino and the atmospheric neutrino mass 
difference respectively [30| . In the scenario of hierarchical neutrino masses these results suggest mi ~ 0, m 2 ~ Am so i ar , 
and m 3 ~ Ar7i atmosphcric . For the inverted hierarchy m 3 ~ 0, m 2 ~ Am atmospheric , and m 1 ~ Ara atmosphoric . However, 
if the neutrino masses are degenerate, one finds mi ~ to 2 ~ 7713 3> Am a t moS pheric- These oscillation experiments can 
merely constrain the squared mass differences, Am 2 , they cannot detect the absolute value of neutrino masses which 
'■^ I is of great significance. 

However, cosmological observations can put upper limits of the absolute neutrino mass. For background evolution, 
neutrino masses contribute to the cosmic energy budget so that modify the epoch of matter-radiation equality, angular 
+3 , diameter distance to the last scattering surface (LSS) and other related physical quantities. For the perturbation, 
neutrinos become non-relativistic at late time thus they damp the perturbation within their free streaming scale. This 
suppressed the matter power spectrum by roughly AP/P ~ — 8Q v /£l m [3lj • Dark energy can also affect the evolution 
of background and perturbation, which can mimic the behavior of neutrino to some extent, so that there exists an 
k>( I obvious correlation between equation of state of dark energy and neutrino mass. 

For simplicity, Ref. [H, HI] parameterizes dark energy by constant equation of state and finds it strongly correlates 
with neutrino mass. In this paper, we study this correlation with dynamical dark energy models which are more 
generic with current cosmological observations as well as with future simulated data. Further we have determined the 
parameters of dynamical dark energy and neutrino mass limit simultaneously. 

The rest part of this paper is structured as follows: in the next section we present our fitting method and data we 
use. In section III we give our results. We make conclusion and discussion in the last section. 



II. METHOD AND DATA 

We choose the commonly used parametrization of the EOS of dark energy as [3^ ]: 

V 

w(z)^w + w 1 — — , (1) 
1 + z 

where z is the redshift. Clearly the cosmological constant corresponds to wq = — 1, w\ = thus ACDM model is 
incorporated in this parametrization. For neutrino mass, we use the parameter /„ which is defined as dark matter 
neutrino fraction: 

U=Pu/PD M - Q3AQ5 e y Qch2 , (2) 
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where p v and pdm denote the energy density of neutrino and dark matter respectively, Y,m v is the sum of neutrino 
mass and Q c h 2 is the physical cold dark matter densities relative to critical density. 

The fitting and statistics strategy we adopt is based on the Markov Chain Monte Carlo package CosmoMC [35]] 1 . 
which has been modified to implement dark energy perturbations with EOS getting across —1 [1 91 ] - Our most general 
parameter space is 

p = (u) b , uj c , Q s , t, w , wi, /„, n s , logflO 10 ^]) , (3) 

where Ub = fl^h 2 and lo c — il c h 2 are the physical baryon and cold dark matter densities relative to critical density, Qs 
is the ratio (multiplied by f 00) of the sound horizon to the angular diameter distance at decoupling, r is the optical 
depth, A s is defined as the amplitude of initial power spectrum and n s measures the spectral index. Assuming a flat 
Universe motivated by inflation and basing on the Bayesian analysis, we vary the 9 parameters above and fit to the 
observational data with the MCMC method. We take the weak priors as: r < 0.8, 0.5 < n s < 1 .5, —3 < wo < 3, — 5 < 
w\ < 5 2 , < f u < 0.5 and a cosmic age tophat prior as 10 Gyr< t <20 Gyr. Furthermore, we make use of the HST 
measurement of the Hubble parameter Ho = lOO/i km s _1 Mpc~ 1 [36] by multiplying the likelihood by a Gaussian 
likelihood function centered around h = 0.72 and with a standard deviation a = 0.08. We impose a weak Gaussian 
prior on the baryon density flbh 2 = 0.022 ± 0.002 (la) from Big Bang nucleosynthesis [37j |. 

In our global fitting we compute the total likelihood to be the products of separate likelihoods of CMB, SNIa, LSS 
and Heidelberg-Moscow experiment (HM) [H, [H, H(| • Alternatively defining \ 2 — ~21og£, so 

Xtotal = XCMB + XSNIa + XlSS + XhM ■ (4) 

For CMB data we use the three-year WMAP (WMAP3) Temperature- Temperature (TT) and Temperature- 
Polarization (TE) power spectrum with the routine for computing the likelihood supplied by the WMAP team [ill ]. 
The supernova data we use are the "gold" set of 157 SNIa published by Riess et aim. 0]. We have marginalized over 
the nuisance parameter [42 in the calculation of SN la likelihood. For LSS information, we have used the 3-D matter 
power spectrum of SDSS [4] and 2dFGRS §, Lyman-a forest data (Lya) from SP SS [43| and recent measurement of 
the baryon acoustic oscillation feature in the 2-point correlation function of SDSS [44|. To be conservative but more 
robust, we only use the first 14 bins of the SDSS 3-D matter power spectrum, which are well within the linear regime 
[45l | . For Lya likelihood, we modify the interpolating code 3 to incorporate dynamical dark energy models. For BAO 
likelihood, we use the constraint 4 |44| : 

A = D v (().ih) = 0,469 ± 0.017 , (5) 
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where H (z) is the Hubble parameter, c is the speed of light and Dm(z) is the comoving angular diameter distance at 
a specific redshift z. Thus 



2 _ (A - 0.469) 2 

Xbao - 



0.017 2 



(7) 



Moreover, the Heidelberg-Moscow experiment uses the half time of 0v2(3 decay to constrain the effective Majorana 
mass and this translates to the constraints of sum of neutrino mass under some assumptions [ttT ] : 



T 1 m v - 1.8 ± 0.6eU (2a) 



(») 



1 Available at http://cosmologist.info/cosmomc/ 

2 We set the prior of wo and wi broad enough to ensure the EoS can evolve in the whole parameter space. 

3 Available at http://www.cita.utoronto.ca~pmcdonal/LyaF/public.lyafchisq.tar.gz 

4 In this work we mainly focus on the correlation between dark energy parameters and the neutrino mass rather than the absolute value of 
the neutrino mass. What we are interested in is the possible effect of BAO measurement on this correlation. Since the BAO measurement 
seems not very consistent with the other observations, such as CMB and SN la, on the constraints of dark energy parameters, we don't 
use the full power spectrum analysis of BAO in our calculations for simplicity. Similarly the author of Ref . |46| also use the parameter 
A to study the neutrino mass limit. 
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Given the Heidelberg-Moscow experiment is controversial for the time being and seems not consistent with other 
observational data, we just make a tentative fit choosing the HM prior to illustrate the effect on the correlation 
between the dark energy parameters and the neutrino mass when the total neutrino mass is very large. 

To get robust conclusion, we have also used the simulated data from future observations of PLANCK, LAMOST 
and SNAP. The fiducial model we choose is the best-fit from the WMAP3 dataset [l|: 6 = 0.04, fl c = 0.20, ho = 0.73, 
r = 0.088, n s — 0.95 and A s = 0.68. And we assume the sum of neutrino mass — and the equation of state 

of dark energy Wq = — 1, w% = 0. 

For the CMB simulation we consider a simple full-sky (f s k y = 1) simulation at Planck- like sensitivity and ignore 
the lensing effect and the tensor information. We neglect foregrounds and assume the isotropic noise with variance 
Nf T = N EE /2 = 3 x 10 -4 /j,K 2 (Pessimistic Planck-like sensitivity) and a symmetric Gaussian beam of 7 arcminutes 
full- width half-maximum (FWHM) (48[. We use the simulated Cj up to I — 2500 for temperature and I — 1500 for 
polarization. The effective x 2 is: 

2 o, r ^/ n , ,J, (cJ T Cf E -{Cj E f\ Cf T C EE + Cf T C EE ~2Cj E Cj E 1 

where C^ Y denote theoretical power spectra and C* Y denote the power spectra from the simulated data. The 
likelihood has been normalized with respect to the maximum likelihood, where Y = C^ Y [49|, [5(| ■ 

For the future LSS survey we consider LAMOST project. The Large Sky Area Multi-Object Fiber Spectroscopic 
Telescope (LAMOST) project as one of the National Major Scientific Projects undertaken by the Chinese Academy 
of Science, aims to measure ~ 10 7 galaxies with mean redshift z ~ 0.2 [HTj] . In the measurements of large scale matter 
power spectrum of galaxies there are generally two statistical errors: sample variance and shot noise. The uncertainty 
due to statistical effects, averaged over a radial bin Afc in Fourier space, is [52| 

The initial factor of 2 is due to the real property of the density field, V is the survey volume and n is the mean galaxy 
density. In our simulations for simplicity and to be conservative, we use only the linear matter power spectrum up to 
k ~ 0.1 h Mpc- 1 . 

The projected satellite SNAP 5 (Supernova / Acceleration Probe)would be a space based telescope with a one square 
degree field of view with 1 billion pixels. It aims to increase the discovery rate for SNIa to about 2J300 per year [53|. 
The simulated SN la data distribution is taken from Refs.[54],[55|]. As for the error, we follow Ref.[54[ which takes the 
magnitude dispersion 0.15 and the systematic error <j sys (z) — 0.02 x z/1.7, and the whole error for each data is 



= l/< a (*) + ^T- . (11) 



where rij is the number of supernova in the i'th redshift bin. 

For each regular calculation, we run 6 independent chains comprising of 150,000-300,000 chain elements and spend 
thousands of CPU hours to calculate on a cluster. The average acceptance rate is about 40%. And for the convergence 
test typically we get the chains satisfy the Gelman and Rubin [56j criteria where R-l<0.1. 

Despite our ignorance of the nature of dark energy, it is more natural to consider the DE fluctuation whether DE 
is regarded as scalar field or fluid rather than simply switching it off. The conservation law of energy reads: 

T> a, . ll = , (12) 

where T MI/ is the energy-momentum tensor of dark energy and ";" denotes the covariant differentiation. Working in 
the conformal Newtonian gauge, Equation (fT2"|) leads to the perturbation equations of dark energy as follows [57l ]: 

8 = -(1 + w)(9- 3$) - 3H(8p/8p - w)S , (13) 

= -H(l - 3w)6 - -^6 + k 2 (^^S + *) . (14) 
1 + w 1 + w 



5 SNAP is one of the several candidates emission concepts for the Joint Dark Energy Mission(JDEM). 
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TABLE 1. Median values and la constrains on dark energy parameters and sum of neutrino mass for models discussed in the 
text. From left to right, we use WMAP three year data (WMAP3), WMAP3+Riess 157 "Gold" 
sample+2dF+SDSS(14bands)(LSS), WMAP3+LSS+Lyman-a forest(Lya), WMAP3+LSS+Lyman-a+Baryon Acoustic 
Oscillation (BAO) information, WMAP3+LSS+Lya+BAO+Heidelberg-Moscow experiment data respectively. The right 
most column we use the simulated data for future SNAP+PLANCK+LAMOST. For one-tailed distributed parameter such as 

neutrino mass, we quote the 95% C.L. limit. 
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-0.067 
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m„+Dynamical Dark Energy without Dark Energy Perturbation 
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< 0.433 


0.899±g;IH 
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-0.131 


-0.236 
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For the models where the EOS doesn't cross -1, the above equation JT3]), (fT4"l) is well defined. For the crossing models, 
the perturbation equation (flU]) , ([14")) is seemingly divergent. However basing on the realistic two-field-quintom model 
as well as the single field case with a high derivative term the perturbation of DE is shown to be continuous 
when the EOS gets across -1, thus we introduce a small positive parameter e to divide the full range of the allowed 
value of the EOS w into three parts: 1) w > — 1 + e; 2) — 1 + e > w > — 1 — e; and 3) w < — 1 — e. 

For the regions 1) and 3) the perturbation is well defined by solving Eqs. fTS")) . (fl4")) as shown above. For the 
case 2), the perturbation of energy density 5 and divergence of velocity, 9, and the time derivatives of 5 and 9 are 
finite and continuous for the realistic quintom dark energy models. However for the perturbations with the above 
parameterizations clearly there exists some divergence. To eliminate the divergence typically one needs to base on 
the multi-component DE models which result in the non-practical parameter-doubling. A simple way out is to match 
the perturbation in region 2) to the regions 1) and 3) at the boundary and set [ljj, |2Cj, |2l| 

5 = , = 0. (15) 

We have numerically checked the error in the range \Aw = e| < 10 -5 and found it less than 0.001% to the exact 
multi-field quintom model. Therefore our matching strategy is a perfect approximation to calculate the perturbation 
consistently for crossing models. For more details of this method we refer the readers to our previous companion 
papers [11 HjJ, HI . ' 

III. RESULTS 

In this section we present our global fitting results of dark energy parameters wq, w\ and neutrino mass limit. We 
pay particular attention to the correlation between w\ , which delineates the time evolving of dark energy, and neutrino 
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mass. Our main result was summarized in Table I. We choose the general model of 9 free parameters as illustrated 
in Eq. §3§ to study the correlation between dynamical dark energy and neutrino mass and discuss the model with DE 
with constant equation of state (WCDM) and ACDM model for comparison. Here we only list the parameters of dark 
energy and neutrino mass. Since dark energy perturbation is crucial in determination of cosmological parameters 
especially for dark energy parameters pi [l^. I20I |22| . we include the full dark energy perturbation in our global fitting 
as well as showing the result by incorrectly switching off dark energy perturbations to reinforce this key point. We 
use different data combination discussed in previous section. The cov(X,Y) is the correlation coefficient of samples 
which quantifies the correlation between the two parameters X and Y defined as [581]: 



cov(X,Y) 



Er=i(^-*i)V£r=i( y *-^) 



(16) 



where the bar denotes the mean value of parameters. If cov(X,Y)>0 means X,Y is positively correlated and vice 
versa. 
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FIG. 1: 95% C.L. contour plots of the sum of neutrino mass Em„ with dark energy parameters wo, wi as well as as 
with different combination of datasets. From outside in , we use WMAP3 only (black), WMAP3+RIESS+SDSS+2dF 
(red), WMAP3+RIESS+SDSS+2dF+Lya (cyan), WMAP3+RIESS+SDSS+2dF+Lya+BAO (purple) and simulated data of 
SNAP+PLANCK+LAMOST (innermost blue shaded). The two lines in the u>o, wi plane distinguish the dark energy models 
(see text) and their intersecting point denotes the ACDM model. In the fitting we include the full dark energy perturbations 
of dynamical dark energy models. 



From the comparison of results with/ without dark energy perturbation, we find dark energy perturbation modifies 
the mean value of all the parameters as well as enlarging the corresponding errors which is consistence with previous 
analysis in the literature [E IS [20, S • We will merely discuss the results with dark energy perturbation hereafter. 
We have seen that WMAP3 only put weak constraints on wo, wi and neutrino mass since CMB data is sensitive to 
the effective equation of state of dark energy defined as [59| : 

_ JdaQ(a)w{a) 
Weff = Jdan(a) ' 

thus it's hard to constrain wq, w\ and neutrino mass by CMB data alone. Adding SN la and LSS data we find the 
mean value of each parameter gets modified and the error bars shrink a lot. However it is noteworthy that adding 
BAO data enlarges the error of wq, w\ while tightens the neutrino mass limit to a great extent. If we take the HM 
prior in Eq.®, we find the error of wq, w\ become greater while the 1-D posterior distribution of neutrino mass turns 
to be two-tailed. We have also done the global fitting using the simulated data of future observation and found the 
errors shrink greatly. 
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FIG. 3: Evolution of the equation of state of dark energy constrained by different combination of datasets as illustrated in 
the plots. The blue solid curve and the gray shaded bands show the mean value and 1-a error respectively. Dark energy 
perturbation is fully included. 
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In Fig. ([I]) we show the 95% C.L. 2-D contours of wq, w%, neutrino mass and as. From outside in, we add data as 
illustrated in the figure caption. Using simulated data we obtain the inner most shaded regions. The contours shrink 
significantly while interesting correlations remain. The results including HM prior are illustrated in Fig.([2|). We also 
cut the wo — wi plane into four parts for different dark energy models as we did in [2(| [U • The EOS of Quintom 
A crossing -f from upside down, i.e. the EOS is greater than -f in the early epoch yet negative than -1 currently 
while quintom B crosses -f from the opposite direction. We find that the ACDM model cannot be ruled out by all 
the data combination at 2a level however the quintom-like dynamical dark energy model is mildly favored. The one 
dimensional constraints of EOS by different combination of datasets are illustrated in Fig. We find that the best 
fit models are always allowed crossing the boundary w = — 1 in these four combinations. When adding the dataset 
of BAO measurement or HM measurement the direction of EOS crossing the boundary will be changed. This is due 
to the changed sign of the mean value of W\ . And the HM prior seems favor more negative value of W\ and stronger 
evolution of EOS. 

In Ref. L 32], the author just used the constant EOS to consider the correlation between EOS of dark energy and 
neutrino mass. The critical point of our result is the correlation among wo, w\ and neutrino mass. From the table 
I we find there exists some correlation of (wo, vriy) and especially for w\ and m„. More interestingly, we find using 
future high quality data, the correlation between wq and m„ nearly vanishes while there is a strong and negative 
correlation between the "running" of dark energy EOS w\ and the total neutrino mass Hm v . Fig. ([4]) delineates the 
fitting result of future simulation in more detail and we clearly see the strong correlations between w\ and neutrino 
mass. In Ref. 60], the authors conclude wi is not strongly correlated with neutrino mass basing on Fisher matrix 
analysis. However, we have seen the distribution of w\ and neutrino mass is highly non-Gaussian and Fisher matrix 
technique is not viable in such cases (50j . Furthermore, they didn't consider dark energy perturbation when EOS 
crosses —1. This also give rise to great bias of determining neutrino mass limit if neglecting time evolving of dark 
energy as illustrated in the table I. 




-1.2 -1 -0.8 -1 -0.5 0.5 0.2 0.4 

w„ w_ Em 

1 v 



FIG. 4: Constraints on dark energy parameters wo, wi and the total neutrino mass Em„ from the Monte Carlo method. The 
two dimensional contours stand for 68% and 95% confidence levels . 

The reason of the correlation between £m„ and dark energy is not hard to understand. The main correlation 
between dark energy EOS and neutrino mass stems from the geometric feature of our universe [32]. In addition, 
dynamical dark energy will modify the time evolving potential wells which affects CMB power spectra through the 
late time ISW effect. Dynamical dark energy can leave imprints on CMB, LSS power spectrum and Hubble diagram 
nonetheless these features can be mimicked by cosmic neutrino to some extent [6 1[. This is the source of correlation. 

IV. CONCLUSION 

In this paper we have done the global fitting including parameters expressing time evolving dark energy, neutrino 
mass and other 6 cosmological parameters using different combination of currently data as well as simulated data 
for future observation. We have seen the interesting correlation between neutrino mass and EOS of time evolving 
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dark energy especially for £m„ and w\. The correlation between the "running" of dark energy EOS Wi and the total 
neutrino mass Em„ is of great significance in determining neutrino mass limit and revealing the nature of dark energy 
by futuristic high precision observational data, such as SNAP, PLANCK 6 and so forth. 
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